Neutrino Anisotropies after Planck 
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We present new constraints on the rest-frame sound speed, 



and the viscosity parameter, 



c^,is, of the Cosmic Neutrino Background from the recent measurements of the Cosmic Microwave 
Background anisotropies provided by the Planck satellite. While broadly consistent with the ex- 
pectations of c^g = Cvis = 1/3 in the standard scenario, the Planck dataset hints for a higher value 
of the viscosity parameter, with c^jg = 0.60 ± 0.18 at 68% c.L, and a lower value of the sound 
speed, with c^g = 0.304 ±0.013 at 68% c.L. We find a correlation between the neutrino parameters 
and the lensing amplitude of the temperature power spectrum Ai^. When the latter parameter 
is allowed to vary, we find a better consistency with the standard model with cJig — 0.51 ± 0.22, 
els = 0.311 ± 0.019 and tIl = 1.08 ± 0.18 at 68% c.L. This result indicates that the anomalous 
large value of Al measured by Planck could be connected to non-standard neutrino properties. 
Including additional datasets from Baryon Acoustic Oscillation surveys and the Hubble Space Tele- 



scope constraint on the Hubble constant, we obtain c: 
Al = 1.15 ±0.17 at 68% c.L. 

PACS numbers: 98.80.Es, 98.80.Jk, 95.30.Sf 
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0.40 ± 0.19, cts = 0.319 ± 0.019, and 



I. INTRODUCTION 

The recent measurements of the Cosmic Microwave 
Background (CMB, hereafter) anisotropies provided by 
the Planck satellite [IHS] are in excellent agreement with 
expectations of the standard ACDM cosmological model 
and present the tightest ever constraints on its parame- 
ters. 

These new observations open the opportunity to fur- 
ther test some of the assumptions of the ACDM model 
and to possibly identify the presence of new physics. 

Following this line of investigation and previous analy- 
ses (see e.g. [^^), in this paper we test some properties 
of the Cosmic Neutrino Background (CNB, hereafter). It 
is well known that, apart from CMB photons, the stan- 
dard cosmological model predicts the existence of a CNB 
of energy density (when neutrinos are relativistic) of 
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where p^ is the photon energy density for T^ = 2.72b K 
and TVeff is the number of relativistic degrees of freedom. 

In the past years, cosmological data increasingly con- 
strained the value of iVcff, ruling out the possibility of 
Ncff = at high significance (e.g. see [3] and references 
therein). More recently, the Planck experiment [3] re- 
ported the bound N^s = 3.51 ± 0.39 at 68% c.L, consis- 
tent in between two standard deviations with the stan- 
dard value Ncs = 3.046 and providing evidence for the 
neutrino background at the level of about nine standard 
deviations. 
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However, the slightly higher value for N^s suggests 
that new physics can be present in the neutrino sector. 
This new physics could also be connected with the lensing 
amplitude of the temperature power spectrum. Let us re- 
mind here that gravitational lensing acts on the CMB by 
deflecting the photon path by a quantity defined by the 
gradient of the lensing potential (p (n) , integrated along 
the line of sight n. Lensing also affects the power spec- 
trum by smoothing the acoustic peaks and this effect 
can be parameterized by introducing the lensing ampli- 
tude parameter Al, as defined in [9], which performs a 
rescaling of the lensing potential 
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where Cf is the power spectrum of the lensing field. In- 
terestingly, the Planck data suggest an anomalous value 
for the lensing amplitude of Al — 1.22^q}3 at 68% c.L, 
i.e. higher respect to the expected value of Al = 1 at 
about two standard deviations. 

Allowing a simultaneous variation in N^g and Al, we 
recently found N^s = 3.71±0.40 and Al = 1.25±0.13 at 
68% c.L in [5], suggesting the presence of some anomalies 
at higher significance. 

In this paper, we continue our search for anomalies 
but considering a different modification to the CNB. In- 
deed, instead of varying the neutrino effective number, 
which we fix at N^s — 3.046, we modify the CNB clus- 
tering properties as first proposed in [TIT . Following [TT] , 
the CNB can be modelled as a Generalized Dark Matter 
(GDM) component with a set of equations, describing 
the evolution of perturbations, given by ffTl lTTlIT^ ): 



q^ = kclg [6„ + 3-^] q^ - -kiT^ 
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When exploring extended cosmological scenarios, we 
vary the effective sound speed c^fj , the viscosity parame- 
ter c^jg and the lensing amphtude parameter Ai^. Firstly, 



where c^g is the sound speed in the CNB rest frame, 
describing pressure fluctuations respect to density per- 
turbations, and cjjg is the "viscosity" parameter which 
parameterizes the anisotropic stress. For standard neu- 
trinos, we have c^g = c^jg = 1/3. Constraints on these 
parameters have been set by several authors (see e.g. 
m [71 [13]), since the observation of deviations from the 
standard values could hint for non-standard physics. Us- 
ing cosmological data previous to Planck, Smith et al. 
[6] found that, assuming iVcff — 3.046, the case c^g = 
c^ig = 1/3 is ruled out at the level of two standard devi- 
ations. More recently, Archidiacono et al. [7j found that 
current cosmological data from the South Pole Telescope 
SPT [14l[T5] exclude the standard value of cl^^ = 1/3 at 
2(7 level, pointing towards a lower value. 

It is therefore extremely timely to bound the values 
of the neutrino perturbation parameters using the recent 
Planck data. Here we provide those new constraints but 
also considering the possible degeneracies between c^g, 
cjjg and the temperature power spectrum lensing ampli- 
tude Al. 

The paper is organized as follows: the next section is 
devoted to the description of the analysis method; our 
results are summarized in Section III; conclusions are 
drawn in Section IV. 



II. DATA ANALYSIS METHOD 

We sample a four-dimensional set of standard cosmo- 
logical parameters, imposing flat priors: the baryon and 
cold dark matter densities il,b and fie, the angular size 
of the sound horizon at decoupling 0, and the optical 
depth to rcionization r. As inflationary parameters we 
consider the scalar spectral index ns and the overall nor- 
malization of the spectrum As at fc = 0.05 Mpc""'^. We 
consider purely adiabatic initial conditions and we im- 
pose spatial flatness. As far as neutrino components are 
concerned, we assume the Planck collaboration |3] base- 
line model: the total number of relativistic degrees of 
freedom A^eff = 3.046, with a single massive eigenstate 
with TTiy = 0.06 eV. We checked that assuming mass- 
less neutrinos does not affect our results. The Helium 
abundance Yp is also varied but assuming a Big Bang 
Nucleosynthesis (BBN) consistency (given Ncs and ilb, 
1^ is a determined function of them) and is therefore not 
treated as an extra free parameter. 



we combine them in pair (c^g - 



^visi ^cff 



Ah, Cyjs — Ah), 



fixing the third parameter at its standard value {Ah — 1, 
•^vis = 1/3) clff = 1/3). Finally, we combine the three pa- 
rameters all together. We adopt the following flat priors: 
0<cli„clg< 1 andO< Al <4. 

Concerning the datasets, we use the high-£ Planck 
temperature power spectrum (CajiiSpec, 50 < ^ < 2500 
J2J) and the \ow-£ Planck temperature power spectrum 
(commander, 2 < i" < 49 [5]), in combination with the 
WMAP low-i likelihood for polarization [l6] at f < 23. 
We refer to this combination as PlanckW P . 

We also consider the inclusion of additional datasets. 
Referring to the latest measurements of the Hubble Space 
Telescope (HST) [T7] , we assume a gaussian prior on the 
Hubble constant Ho = 73.8 ±2.4 km s"^ Mpc"^ We also 
include measurements of Baryon Acoustic Oscillations 
(BAO) from galaxy surveys. Here, we follow the same 
approach reported in Ref. [3j combining three datasets: 
SDSS-DR7 ,18J, SDSS-DR9 [S] and 6dF Galaxy Survey 
|20) . We refer to this combination as PlanckEX. 

Our analysis method is based on the publicly avail- 
able Monte Carlo Markov Chain package cosmomc [5T] 
(version released in March 2013) using the Gelman and 
Rubin statistic as convergence diagnostic. 



III. RESULTS 

We show our results in Tab[l] in the form of the 68% 
confidence level, i.e., the interval containing 68% of the 
total posterior probability centered on the mean. 

As we can see from Tab|Tj FigJT] and Fig[2j allowing 
both Cgfj and cjjg to vary produces posterior values in 
disagreement with the standard model. In particular, 
the constraints on c^jg point towards a value larger than 
1/3, being 4;^ = 0.60 ± 0.18 at 68% c.L, in tension with 
the standard value at about 1.5 standard deviations. On 
the other hand, c^g assumes in our analysis a value lower 
than the expected one: c^g = 0.304 ± 0.013 at 68% c.L, 
ruling out the standard value at more than 95% c.L. The 
addition of BAO and HST datasets only results in an 
improvement of the constraining power. 

The situation is significantly different when the lens- 
ing amplitude parameter is allowed to vary (see Tab[l| 
Figjl] Fig[2] and Fig[|). In the c^j^ + Al case, the 
standard value of the viscosity parameter is recovered 
(cJig = 0.35±0.12); similarly, in the c^^-I-Al case, c^g is in 
agreement with the expected value (c^jj = 0.321 ±0.014). 
However, though the Al = 1 case is now in more agree- 
ment with the data, it is still in disagreement with 
the standard value at more than la level (respectively, 
Al = 1.20 ± 0.12 and Al = 1.16 ± 0.13), showing a de- 
generacy between Al and the clustering parameters c^^g 
and c^jj (see Figjs] and FigJ6| . Also for these choices of 
parameters, the addition of BAO and HST allows to get 



TABLE I. Comparison between 
posterior means for the cosmolog 



extended cosmological models and the standard ACDM for different datasets. Listed are 
iical parameters from the indicated datasets (errors refer to 68% credible intervals). 



Dataset 


Parameter 


ACDM 


±Cvis ± Ccft 


±C?ff ± Al 


±Cvis ± Al 


±C?ff ± cJis + A-L 




100 Qbh^ 


2.206 ±0.028 


2.118 ±0.047 


2.219 ±0.045 


2.236 ±0.053 


2.162 ±0.095 




n,h^ 


0.1199 ±0.0027 


0.1157 ±0.0038 


0.1177 ±0.0032 


0.1170 ±0.0034 


0.1159 ±0.0036 




100 6* 


1.0413 ± 0.0006 


1.0412 ±0.0014 


1.0428 ±0.0012 


1.0421 ±0.0019 


1.0420 ±0.0020 




log[10^''^s] 


3.089 ±0.025 


3.173 ±0.052 


3.086 ±0.028 


3.08 ±0.05 


3.141 ±0.078 




T 


0.090 ±0.013 


0.089 ±0.013 


0.088 ±0.013 


0.087 ±0.013 


0.089 ±0.014 


PlanckWP 


ns 


0.9606 ± 0.0073 


0.998 ±0.018 


0.9732 ± 0.0099 


0.970 ±0.014 


0.989 ±0.023 




A^ 


= 1 


= 1 


1.16 ±0.13 


1.20 ±0.12 


1.08 ±0.18 




Wis 


= 0.33 


0.60 ±0.18 


= 0.33 


0.35 ±0.12 


0.51 ±0.22 




cl^ 


= 0.33 


0.304 ±0.013 


0.321 ±0.014 


= 0.33 


0.311 ±0.019 




Ho ("' 


67.3 ±1.2 


68.0 ±1.3 


68.7 ±1.5 


68.9 ±1.5 


68.6 ±1.7 




100 Qbh'^ 


2.236 ±0.026 


2.142 ±0.048 


2.228 ±0.037 


2.258 ± 0.043 


2.209 ±0.091 




Q,h^ 


0.1180 ±0.0016 


0.1160 ±0.0025 


0.1171 ±0.0018 


0.1173 ±0.0023 


0.1167 ±0.0024 




100 6* 


1.0417 ±0.0006 


1.0417 ±0.0016 


1.0427 ±0.0011 


1.0429 ±0.0018 


1.0426 ±0.0018 




log[10^«As] 


3.111 ±0.032 


3.158 ±0.050 


3.084 ±0.028 


3.065 ± 0.043 


3.104 ±0.075 




T 


0.101 ±0.015 


0.091 ±0.013 


0.088 ±0.013 


0.087 ±0.013 


0.088 ±0.013 


PlanckEX 


ns 


0.9621 ±0.0061 


0.993 ±0.016 


0.9737 ±0.0080 


0.967 ±0.012 


0.978 ±0.023 




Al 


= 1 


= 1 


1.18 ±0.12 


1.23 ±0.11 


1.15±0.17 




Wis 


= 0.33 


0.53 ±0.16 


= 0.33 


0.302 ±0.097 


0.40 ±0.19 




Ceff 


= 0.33 


0.306 ±0.013 


0.322 ±0.013 


= 0.33 


0.319 ±0.019 




^0 ("' 


68.31 ±0.73 


68.29 ±0.74 


68.93 ±0.80 


69.16 ±0.87 


68.88 ±0.99 



km s~i Mpc~i 



tighter constraints as welL 

Finally, when all the three parameters are allowed to 
vary, their posteriors are in good agreement (within la 
c.l.) with the standard model: we find c^g = 0.311 ± 
0.019, c^i^ = 0.51 ± 0.22 and Al = 1.08 ± 0.18. In this 
case, the addition of BAO and HST datasets slightly 
affects the peak location in the posterior distributions 
(see Tabjl) Figjl] Fig[2]and Fig[3|. In particular, it pro- 
duces a 0.5ct displacement of the c^jg mean towards a 
smaller (therefore more compatible with the expected) 
value. Conversely, it results in a 0.3 a shift of the Al 
mean towards a higher value (less compatible with the 
expected, but more compatible with the Planck collabo- 
ration finding). On the other hand, the magnitude of the 
displacement induced on the c^g distribution (towards 
the expected value) is minimal (< 0.3 cr). 



have found that the Planck dataset hints for anomalous 
values for these parameters with c^j^ = 0.60±0.18 at 68% 



c.l. and 



'off 



0.304 ± 0.013 at 68% c.l. 



We found a correlation between the neutrino parame- 
ters and the lensing amplitude of the temperature power 
spectrum A^. When this parameter is allowed to vary 
we found a better consistency with the standard model 



with 



0.51 ± 0.22, c; 



eff 



0.311 ± 0.019, and 



Al = 1.08 ± 0.18 at 68% c.l.. This resuh indicates that 
the anomalous large value of Al measured by Planck 
could be connected to non-standard neutrino properties. 
Including additional datasets from Baryon Acoustic 
Oscillation surveys and the Hubble Space Telescope con- 
straint on the Hubble constant we obtain c^jg = 0.40 ± 
0.19, clff = 0.319 ± 0.019, and Al = 1.15 ± 0.17 at 68% 
C.L. 



IV. CONCLUSIONS 
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FIG. 1. One-dimensional posterior probabilities of the parameter c^g for the indicated models for PlanckWP (left) and 
PlanckEX (right) datasets. The vertical dashed line indicates the expected value in the standard model. 




FIG. 2. One-dimensional posterior probabilities of the parameter cJig for the indicated models for PlanckW P (left) and 
PlanckEX (right) datasets. The vertical dashed line indicates the expected value in the standard model. 
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FIG. 3. One-dimensional posterior probabilities of the parameter A^ for the indicated models for PlanckW P (left) and 
PlanckEX (right) datasets. The vertical dashed line indicates the expected value in the standard model. 
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FIG. 4. Two-dimensional posterior probabilities in the cJis — Cpg plane for the indicated datasets and models. The dashed lines 
indicate the expected values in the standard model. 
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FIG. 5. Two-dimensional posterior probabilities in the Ccs — j4l plane for the indicated datasets and models. The dashed lines 
indicate the expected values in the standard model. 
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